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Abstract 

We show that the tails of the single-point velocity probability distribution 
function (PDF) are generally non-Gaussian in developed turbulence. By us- 
ing instanton formalism for the Navier-Stokes equation, we establish the re- 
lation between the PDF tails of the velocity and those of the external forc- 
ing. In particular, we show that a Gaussian random force having correlation 
scale L and correlation time r produces velocity PDF tails lnV{v) oc —v'^ at 
V ^ v-rras^L/r. For a short-correlated forcing when r <C L/vrms there is an 

intermediate asymptotics lnV{v) oc —v"^ at L/t ^ v ^ ?^rms- 
PACS numbers 47.10.-Fg, 47.27.-i, 05.40.-Fj 
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Early experimental data on skewness and flatness of the velocity fleld prompted one 
to believe that the single-point velocity PDF in developed turbulence is generally close to 
Gaussian |1|J^. A possible reasoning may be that large-scale motions (that give the main 
contribution into velocity statistics at a point) are connected to a random external forcing 
/ then velocity v{t) = /q f{t')dt' has to be Gaussian if t is larger then the correlation time r 
of the forcing, irrespective of the statistics of /. That would be the case if the force was the 
only agent affecting velocity. However, there are also nonlinearity (leading to instability and 
break-up of large-scale motions) and viscosity (that dissipates small-scale modes appeared as 
a result of the instability). Let us flrst explain the simple physics involved and formulate the 
predictions following from physical arguments, then we develop the formalism which gives 
the predicted PDF tails. Qualitatively, one may describe the interplay between external 
force and nonlinearity in the following way. Force / pumps velocity v ^ ft until the time 
t* ~ L/f when nonlinearity restricts the growth. The relation between velocity and forcing 
can thus be suggested as follows: f ^ ~ f'L. Therefore, velocity's PDF can be obtained by 
substituting / ~ v'^/L into force's PDF Vf: V{v) ~ Vf{v^/L). Those arguments presume 
that is less than the correlation time r of forcing. If opposite is the case ^ t then the 
law of velocity growth is different f ^ ~ /"^tr, so that the velocity increases up to ~ f'^Lr; 
the short-correlated pumping is effectively Gaussian Vf ~ exp (-/2) and the velocity's PDF 
is V{v) ~ exp[— (f/frms)^]- We see that velocity PDF is expected to be dependent on the 
statistics of the force. 

Actual mechanism of restriction for the Navier-Stokes equation (instability of a large- 
scale flow leading to a cascade that provides for a viscous dissipation) is irrelevant for the 
above arguments. What matters is that we deal with the system of the hydrodynamic type 
so that nonlinear time t^, can be estimated as L/v. For example, the same argument goes 
for Burgers equation where t^, is a breaking time and viscous dissipation at a shock prevents 
further growth 0. It is interesting that viscosity does not enter above estimates yet it's 
existence is crucial for the whole picture to be valid. 

Let us stress that the above arguments can be only applied to rare events with velocity 
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and force being much larger than their root-mean-square values when the influence of back- 
ground fluctuations can be neglected. The above predictions are thus made for PDF tails. 
For a a nonlinear dissipative system, it is generally difficult to relate an output statistics 
to the statistics of the input (be it initial conditions or external force). Our aim here is 
to show that it is possible, nevertheless, to relate the probabilities of rare fluctuations that 
is to relate the tails of the PDFs of the force and the field that is forced respectively. A 
rigorous way to describe rare fluctuations is the instanton formalism recently developed for 
turbulence P] and employed for obtaining PDF tails in different problems 0-^]. 

The main idea of the method is that the tails are described by saddle-point configurations 
in the path integral for the correlation functions of the turbulent variable (say, velocity v). 
We call the configuration instanton because of a finite lifetime. One may call it also optimal 
fluctuation since it corresponds to the extremum of the probability. 

We start with the Navier-Stokes equation 



where / is a random force (per unit mass) pumping the energy into the system and u is 
the viscosity coefficient. Incompressibility is assumed so that diw = div/ = 0. The 
field P in (P is the pressure divided by the mass density p. Velocity correlation functions 
can be presented as path integrals which form is determined by the statistics of pumping. 
Let us first consider a Gaussian forcing with the correlation function (/^(t, r)/^(t', r')) = 
Sa/3(t — t', r — r') which is assumed to decay on the scale r as a function of the first argument 
and on the scale L as a function of the second one. Then moments of the velocity can be 
written as path integrals: 



(1) 




(2) 



where p is an auxiliary field and the effective action has the following form 0. 





(3) 
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The independent fields P and Q play the role of Lagrange multipliers enforcing the incom- 
pressibility conditions on the velocity and the response field: VaVa = ^aPa = 0. 

The tails of the velocity PDF are determined by high moments with n ^ 1 which can be 
found by applying the saddle-point method to the integral (0): (f^") = u"'{0,0) exp(Xextr)- 
The configuration u(t, r), p(t, r), P(t, r) and Q{t, r) corresponding to a saddle point is our 
instanton. The extremum conditions for the argument of the exponent in @) determining 
the instanton give two dynamical equations 

= -i j dt' dr' Eapit -t',r- r')pf3{t\ r') , (4) 

= 2in^(t)5(r)M„/^^^ (5) 
and two incompressibility constraints 

V'P = -V^iu^V^Uo) . (6) 

V^Q = Va{P(NaUp - UpVpPa) ■ (7) 

As it was explained in f^, the system (§-0) is to be solved at negative time because it 
describes the prehistory that leads to a given measured value of the velocity at t = 0. Since 
p = at t > then the right-hand side of (||) directly produces the value of p at t = —0: 

'^"/^ W - ^ - ^ j T^TTTTTV- (8) 



Pa(0, r) = —2in 



u2(0,0) 

The equations (||,^ are to be solved at — oo < t < with the constraints (HJ^) under the 
boundary conditions and m— >Oast— oo. Instanton approach thus reduces the statis- 
tical problem (finding PDF tails) to the dynamical problem of finding a particular solution of 
the deterministic equations Note that the instanton equations are deterministic they 

present already substantial simplification with respect to the initial Navier-Stokes equation 
with a random force (which is not suprising because they describe only part of velocity 
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statistics, namely the tails of the single-point PDF). Still, the system (^,^ is quite com- 
plicated and it's complete solution (which will provide an important information about the 
spatio-temporal domains with high velocity) is still ahead of us. The main difficulty is an 
effective spatial nonlocality of the constraints 

Our goal in this paper is to show that a definite statement about the functional form 
of velocity PDF tails can be obtained from the analysis of the symmetries of the instanton 
equations without actually finding a solution. Let us show how the symmetry analysis 
gives n-dependence of (t;^") which determines the functional form of PDF tails. What is 
important to state is that both fields u and p decay at moving backwards in time during 
some characteristic time which we call the lifetime of the instanton. Actually it is the 
same time which we qualitatively discussed above at treating rare events with strong forcing 
since the instanton presents just the space-time picture of those typical events contributing 
to (m^"). Symmetry analysis depends on whether the pumping correlation time r is larger 
or smaller than the instanton lifetime t^,. Let us ffist consider the case t ^ t^,, that makes it 
possible to consider the pumping correlation function S as time independent. In this case, 
the parameter n can be excluded from (^-|^ by the rescaling transformation 

t^X'H, u^Xu, P^X^P, Q^X^Q, 

u-^Xu, X^p , = n . (9) 

That gives a general n-dependence of the velocity 

u = n^/^cp{iyyn) , (10) 

where dimensionless function (p is expected to go to some constant when it's argument goes 
to zero. This is equivalent to the physically plausible assumption that the high velocity 
moments are viscosity independent. Under such an assumption the n-dependence of the 
instanton solution u oc n^^^ gives the following n-dependence of the moment (f oc n"/^ 
which corresponds to the PDF tail 

\nV{v^)oc-v\ (11) 
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Note that the integral term in the action lextr oc n i.e. the factor exp{Iextr) gives only sub- 
leading contribution into (p!l]). If L/r < t>j.ms (where frms is the typical value of the velocity 
fluctuations) then the asymptotics ( pH]) is realized at f ^ frms and InP ~ —{v/vj-jas)^- In the 
opposite limit of a short-correlated pumping with L/r Vrms, an intermediate asymptotics 
exists where the pumping correlation function S can be treated as delta-correlated in time: 
S()f:,r) = 5(t)x{r). Then the only changes in (P) are 

p^X^p, Q^X^Q, X^ = n. (12) 

That leads to the law (cf. 0) 

\nV{v')r~.-{v/v„^,y, (13) 

which is valid at L/r ^ f ^ Vy-ms- For larger v, the asymptotics (|TT]) is realized. We thus 
see that for Gaussian pumping the velocity PDF always decreases faster than Gaussian. 

Note that the above consideration can be readily extended for the consideration of 
velocity differences w{r) = \v{r) — f(0)| at any r. There is an intermediate region 
'U^rms(^) ^ w{r) ^ frms wherc V{w) is not described by the direct instanton formalism. The 
same formulas ( [TI|JTB| ) as well as below ( PPD describe the remote tails of V{w) at w v^-^s, 
which are probably unaccessible experimentally at present. 

Contrary, faster-than-Gaussian decay of a single-point velocity PDF was recently ob- 
served both in experiments M and in numerical simulations |[9|-[TT|, the qualitative reason 



for that (short life time of strong fluctuations) was discussed in |]TT|. Unfortunately, except 
the qualitative statement on faster-than Gaussian decay there were no definite data in [p|-[TT 
to allow for a quantitative comparison with (pJ] , p!3D , this remains for the future work. 



The transformations (p|,p!2D shows that the lifetime of the instanton is short for large 
n: oc n"^/^ for the fast pumping r <^ and t^, oc n~^^'^ for the slow pumping t ^ t^,. 
Therefore at large enough n we always deal with a slow pumping. In other words, distant 
tails of the velocity PDF are determined by a simultaneous statistics of the pumping (which 
is not necessary Gaussian) rather than by integral over time as one would naively expect. 
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As we shall show now, the same property (r ^ t^) is correct for a non-Gaussian statistics 
of the pumping. Indeed, the above procedure can be readily generalized for an arbitrary 
pumping statistics when instead of (0) one has 

(|^;|2") = JvfVpVvVPVQVfexp(ti + 2nln\v\) . (14) 

Here, we substituted the simultaneous PDF of pumping Vf{f) since we treat the slow 
pumping and the effective action I figuring in ( p^ for the Navier-Stokes equation has now 
the following form 

J = J dtdr [paidtVa + pVa - I'V'^Va + VqP) 

+ Q'^aVa - Pafa] ■ (15) 

The saddle-point equation for the instanton velocity u is now written as 

dtUa + UfsVf^Ua - uV'^Ua + Va-P = fa , (16) 

with the relation 

J dtpo, = -t6lnVf{f)/6U. (17) 

The equation for p and the incompressibility constraints are the same (^,^,0). 

Let us assume the tail of Ai = —\nVf to be scale invariant with the exponent a: 

M{Xf)=X-M{f), (18) 

Now, we can generalize (|^) 

t^X-h, u^Xu, P^X^P, 

v^Xv, f^X'f, (19) 
p-^X'^^-'^p, Q-^X^^'Q, X2'^ = n, 



keeping the form of the equations ((T^), ([T^) and removing n from the equation (|^). Using 
the law (^) we can find how the characteristic velocity scales with n. As above the answer 



can be expressed in terms of the behavior of the tail of the simultaneous PDF V{v'^) for the 
velocity: 

\nV{v^) (X -v^" . (20) 

The lifetime of our instanton scales as oc n^^l^"^"-^ that is decreases with n for any 
positive a which justifies our consideration. For Gaussian statistics a = 2 and we reproduce 
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Note that the velocity PDF decays always faster than that of the force; in particular, 
Gaussian velocity would correspond to the force PDF decaying exponentially. 

The expression implies homogeneous pumping. Meanwhile, the transformations 



(P,|l^,|T^) do not transform coordinates and can be generalized for spatially ingomogeneous 
pumping statistics. Therefore our conclusions are true also for a more physical case of an 
inhomogeneous pumping, in particular acting on the boundaries of the flow. 

To conclude, let us describe the status of the results obtained: Under the assumptions 
that the solutions of the instanton equations exist and the probability of finding very high 
velocity is viscosity independent we found how that probability is related to the statistics of 
the force. 
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